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Gentlemen: 

The  extension  of  Contract  AF  33(657)- 16678  is  still  operating 
under  a  $250,  000  letter  contract.  Approval  of  the  one  year  extension  in  the 
amount  of  $490,  023  is  in  progress. 

Below  are  set  forth  the  work  progress  against  the  major  tasks 
remaining  tuider  this  contract. 

PUBLICATION  OF  SPECIAL  REPORTS  DEALING  WITH  WORK  COMPLETED 
UNDER  ORIGINAL  CONTRACT 


The  following  special  reports  have  been  published. 

STATEMENT  #2  UNCLASSIFIED 

This  document  is  subject  to  special  export  controls  and  each 
r/uismittal  to  foreign  government  s.  op  foreign  nationals  mty.be 
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LASA  Special  Report  No.  1  -  A  STUDY  OF  THE  RELATIVE 
CAPABILITY  OF  LARGE  AND  SMALL  SEISMIC  ARRAYS 
FOR  EVENT  IDENTIFICATION. 

LASA  Special  Report  No.  2  -  RESEARCH  ON  HIGH 
RESOLUTION  FREQUENCY-WAVENUMBER  SPECTRA. 

LASA  Special  Report  No.  5  -  DETECTION  OF  DISCRETE 
ARRIVALS  IN  THE  MANTLE  P-WAVE  NOISE. 

LASA  Special  Report  No.  8  -  SHORT  PERIOD  SIGNAL 
SIMILARITY. 

LASA  Special  Report  No.  9  -  LONG  PERIOD  SIGNAL 
SIMILARITY. 


The  following  reports  have  been  submitted  to  Captain  Lam  for 
approval  and  are  in  the  final  stages  of  publication. 


LASA  Special  Report  No.  3  -  SUBARRAY  PROCESSING. 

LASA  Special  Report  No.  6  -  ANALYSIS  OF  SUBARRAY 
WAVENUMBER  SPECTRA. 

LASA  Special  Report  No.  12-ANALYSIS  OF  LONG  PERIOD 
NOISE. 

LASA  Special  Report  No.  13-NOISE  COHERENCE  AMONG 
SUBARRAYS. 


The  following  reports  have  not  yet  been  submitted  to  Captain 
Lam  and  are  in  various  stages  of  typing  and  editing. 

LASA  Special  Report  No.  4  -  SPACE  AND  TIME  VARIABILITY 
OF  THE  LASA  NOISE  FIELD. 

LASA  Special  Report  No.  10-EQUALIZATION  STUDIES. 

LASA  Special  Report  No.  11-RESOLUTION  OF  EVENTS. 

LASA  Special  Report  No.  14-WIENER  NON-TIME 
STATIONARY  PROCESSING. 
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CONTINUED  SHORT  PERIOD  NOISE  ANALYSIS 

Three  additional  short  period  noise  tapes  have 
been  received  but  have  not  yet  been  transferred  into  the  Texas 
Instruments'  system. 

The  continued  short  period  noise  analysis  will  con¬ 
sist  of  a  study  of  locally  generated  noise  and  k-line  spectral  analy¬ 
sis  of  a  few  summer  noise  samples. 

At  the  LASA  the  day  noise  is  characteristically 
richer  in  higher  (2  -  3cps)  frequencies.  This  is  almost  certainly 
locally  generated  noise.  A  study  will  be  undertaken  to  describe 
this  locally  generated  noise  as  completely  as  possible.  A  few 
examples  of  this  noise  type  are  already  in  our  short  period  li¬ 
brary  and  these  should  be  sufficient  for  this  study.  This  portion 
of  the  short  period  noise  study  has  not  yet  been  clearly  defined. 
Some  visual  analysis  of  the  noise  samples  will  be  necessary  be¬ 
fore  finalizing  the  study  plans . 

About  three  to  five  summer  noise  samples  will 
be  analyzed  using  interpolated  k-line  spectra.  W^en  considered 
with  the  previous  results  obtained  in  winter  noise  ,  these  should 
give  an  indication  of  seasonal  variations  in  the  major  noise  types. 

LONG  PERIOD  NOISE  ANALYSIS 


The  long  period  noise  analysis  has  been  broken 
down  into  four  tasks. 

•  Stability  of  Noise  Statistics 

Two  signal  extraction  filters  and  one  prediction 
filter  have  been  designed  from  one  noise  sample 
and  applied  to  a  noise  sample  taken  6  hours  later. 


* 


LASA  Special  Report  No.  6  , 

Analysis  of  Subarray  Wavenumber  Spectra 
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Most  of  processing  has  been  completed  and  results 
are  presently  being  analyzed. 

e  Further  Study  of  the  Noise  Below  ■  05  cps 

Previous  analysis  indicated  that  the  noise  below 
.05  cps  was  non-seismic.  One  long  noise  sample 
is  presently  in  house  which  will  be  analyzed  in 
detail  in  an  attempt  to  determine  the  nature  of  the 
low  frequency  noise. 

•  Analysis  of  Summer  Noise  Samples 

The  noise  samples  used  on  the  previous  contract 
covered  the  period  of  October  to  April.  Several 
noise  samples  recorded  during  the  summer  months 
will  be  compared  with  those  obtained  for  the  winter 
samples.  Analysis  will  consist  of  wavenumber 
spectral  estimates,  multiple  coherence  studies, 
and  frequency  power  spectra. 

•  Study  of  the  Composition  of  the  Noise  Field 

Attempts  will  be  made  to  use  the  3-component 
data  to  define  the  propagation  modes  of  the  noise 
as  a  function  of  frequency.  Specific  analysis 
techniques  are  presently  being  formulated. 

SIGNAL  DISSECTION 

Two  of  the  tasks  relating  to  the  signal  dissection 
problem  have  been  defined  and  initiated. 

•  P-Wave  Extraction 


Techniques  for  obtaining  a  "cleaned-up"  esti¬ 
mate  of  the  P-wave  arrival  (i.  e.  ,  a  noise  and 
scattered-energy  free  P-wave)  are  being 
evaluated.  The  adaptive  processing  scheme  is 
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being  tested  and  results  are  being  compared  with 
those  from  time -shift-an*d- sum  techniques.  The 
purpose  of  this  processing  is  to  reduce  noise  and 
other  signal  phases  while  passing  P-wave  phases. 

•  Summation  Signal  Processing 

A  study  of  the  theoretical  ability  of  the  summation 
to  enhance  various  signal  phases  in  the  presence 
of  ambient  noise  and  interfering  signals  will  be 
undertaken.  This  work  will  be  based  entirely  on 
the  summation  wavenumber  response  of  the  array. 
Summation  processing  will  also  be  evaluated  using 
the  same  events  used  to  evaluate  the  multichannel 
filtering  techniques. 

Other  tasks  relating  to  this  problem  have  not  yet 
been  specifically  outlined.  It  is  anticipated  that  the  tasks  will 
be  defined  by  the  end  of  October  and  initiated  immediately  there¬ 
after.  These  tasks  will  involve  extensive  use  of  the  3-component 
data. 

HIGH  RESOLUTION  SPECTRA  OF  SHORT  PERIOD  SIGNALS 
A.  Work  Completed 

LASA  Travel- Time  Anomaly  and  Moveout  Analysis 

This  analysis  consists  of  a  theoretical  analysis 
of  the  moveout  smoothing  spectral  window  and  a  computation 
analysis  of  the  travel-time  residuals  experienced  at  the  Montana 
LASA.  The  theoretical  analysis  utilized  a  plane  wave  signal 
model  and  a  "boxcar'1  smoothing  function  to  study  the  combined 
effect  of  smoothing  and  move  out  on  frequency  domain  proces¬ 
sing  of  seismic  arrays.  Travel-time  residuals  were  computed 
for  113  small  regions  of  the  earth  using  384  recorded  events. 

Conclusions  reached  in  this  study  were: 

•  The  size  of  a  region  which  can  be  corrected  with 

the  same  travel-time  residuals  varies  with 
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azimuth  and  epicentral  distance. 

©  Both  large  travel-time  residuals  and  moveout- 

smoothing  spectral  window  effects  make  the  use 
of  the  LASA  E  and  F  rings  undesirable  when 
computing  high  resolution  wavenumber  spectra. 

The  details  of  this  analysis  are  contained  in  the 
LASA  Special  Scientific  Report  No.  15,  Analysis  of  Travel-Time 
Anomalies  and  Moveout  at  LASA. 

Study  of  Interative  Techniques  for  the  Solution  of 
Frequency  Domain  Filter  Sets 

The  following  algorithms  were  studied  with  the 
objective  of  improving  the  speed  and  accuracy  of  computing  high 
resolution  wavenumber  spectra. 

•  Steepest  descent  procedure 

•  Finite  interation  procedure 

•  Exact  inverse  procedure 

In  general,  the  exact  inverse  procedure  was 
found  to  be  the  most  efficient  procedure.  The  application  of 
this  technique  to  iterative  updating  and  smoothing  filter  sets  was 
developed. 

Details  of  this  study  are  contained  in  LASA  Special 
Scientific  Report  No.  16,  Interative  Techniques  for  the  Solution  of 
Frequency  Domain  Filter  Sets. 

B,  Future  Plans 

Pilot  Study  to  Investigate  High  Resolution  Wave- 
number  Spectra  Operational  Parameters 


This  study  has  the  following  objectives: 
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®  Determine  the  optimum  transform  gate  overlap  for 

locating  single  events  in  ambient  seismic  noise. 

®  Determine  the  optimum  transform  gate  overlap  for 

locating  two  simultaneous  events. 

©  Select  a  set  of  reference  sensors  which  yield  good 

average  location  capability. 

©  Select  a  method  of  correcting  filter  sets  for  travel¬ 

time  residuals  during  the  computation  of  the  wave- 
number  spectra. 

•  Select  a  method  for  automatically  computing  the 

best  estimate  of  an  event  epicenter. 


Study  of  Adaptive  Procedures  for  High  Resolution  Wave  Number  Spectra 


Adaptive  and  iterative  filter  design  procedures  will 
be  evaluated  using  the  pilot  study  data  ensemble.  This  evaluation 
will  compare  the  location  ability  of  single  segment  F  -  K  spectra 
with  those  generated  by  a  recursive  formula  for  the  inverse  of  a 
covariance  matrix  (a  modification  of  the  exact  inverse  proce¬ 
dure).  If  an  adaptive  procedure  produces  an  improvement  in  the 
locatabiiity  of  event  epicenters,  further  studies  will  be  made  to 
optimize  the  procedure;  and  the  final  processing  will  incorporate 
the  procedure. 

Subarray  Processing  and  Event  Detection 

An  isotropic  MCF  with  a  high  velocity  pass  band 
from  0.  6  Hz  to  1.  6  Hz  will  be  designed  for  each  subarray.  This 
processor  will  be  compared  with  partial  summations  of  subarray 
elements. 

The  on-line  LASA  detection  scheme  is  indicated 
to  be  so  effective  that  we  did  not  seek  new  detection  techniques. 
The  LASA  automatic  event  detector  is  a  simple  scheme  where 
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the  current  peak  energy  in  the  band  pass  0.  95  Hz  to  1. 45  Hz  is  com¬ 
pared  with  the  average  energy  in  that  band  pass  one  minute  previous. 
When  this  ratio  is  greater  than  a  specified  threshold  (empirically 
chosen  to  be  5.  o2)  an  event  is  said  to  be  present.  To  reduce  the  false 
detections  due  to  non-teleseismic  disturbances,  one  automatic  event 
detector  is  attached  to  each  subarray  in  the  LASA  E  and  F  rings.  If 
four  of  these  detectors  detect  an  event  in  a  single  20  second  period,  a 
teleseismic  event  is  said  to  be  present.  The  false  alarm  rate  of  this 
detection  scheme  appears  to  be  2-3  events  per  day  of  continuous  opera¬ 
tion. 


Details  of  the  Automatic  Event  Detector  can  be  found  in  a 
paper  by  H.  W.  Briscoe  and  P.L.  Fleck  titled,  "A  Real-Time  Computing 
System  for  LASA"  found  in  the  Proceedings -Spring  Joint  Computer  Con¬ 
ference,  1966.  A  study  of  event  detection  at  LASA  is  presented  in  MIT 
Lincoln  Laboratory  Technical  Note  1966-26  LASA  On-Line  Detection, 
and  Signal- To-Noise  Enhancement. 

Final  Processing 

A  data  ensemble  consisting  of  39  events  currently  in  TI 
library  will  be  processed  in  a  manner  determined  in  the  pilot  study. 

Several  types  of  composite  events  will  be  generated  and  processed  to  test 
the  location  ability  of  the  high  resolution  F-K  spectra  technique  in  a 
complex  signal  environment.  Other  events,  including  several  recorded 
by  the  E3  extended  subarray,  will  also  be  processed  and  the  results 
analyzed.  Single  subarray  data  will  also  be  processed  when  available. 

EVENT  RESOLUTION 

A  new  series  of  LASA  slow  mode  tapes  arrived  in  the  latter 
part  of  August.  Five  of  these  tapes  contained  LP  teleseismic  events  from 
four  locations:  Virgin  Islands,  Argentina,  Sea  of  Okhotsk,  and  Kurile 
Islands.  Magnetic  tapes  in  a  Texas  Instruments'  compatible  format  have 
been  generated  for  three  of  these  four  events,  with  plots  and  tape  dumps 
of  the  data  now  being  studied.  Transcription  of  the  Kurile  Islands  record¬ 
ing  has  not  been  successful;  after  several  attempts  to  reformat  the  SDL 
data  reel  it  has  been  concluded  that  the  source  magnetic  tape  is  faulty. 

Work  to  datu  indicates  that  these  signals  contain  many  spikes 
most  of  these  spikes  should  be  correctable.  It  has  not  yet  been  ascertained 
how  much  of  this  data  will  be  usable. 
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Some  compound  long  period  events  will  be  generated. 

These  will  be  synthesised  using  real  events  compounded  with  various 
scaling  factors  and  time  shifts  to  simulate  desired  relative  magnitudes 
and  epicenter  locations. 

High  resolution  wavenumber  spectra  will  be  run  on  a 
single  and  some  compounded  long  period  events.  These  spectra  should 
be  useful  in  identifying  phases  masked  by  a  strong,  interfering  event. 

Separation  of  signal  waveforms  from  overlapping  signals 
will  be  attempted  using  theoretical  velocity  filters  and  measured  statistics 
The  ability  to  orient  horizontal  components  will  be  used  in  connection  with 
the  above  techniques.  For  example,  when  trying  to  separate  a  small  sur¬ 
face  mode  wave  masked  by  a  large  one,  taking  the  horizontal  components 
transverse  to  the  large  event  should  give  considerable  rejection.  These 
transverse  traces  could  thus  be  used  as  input  to  velocity  filtering  schemes 
to  further  enhance  the  desired  surface  mode  waves. 

FINANCIAL  STATUS 


The  financial  status  of  the  contract  through  August  1967  is 
summarized  on  the  attached  Cost  Planning  and  Appraisal  Chart. 

Yours  very  truly, 

TEXAS  INSTRUMENTS  INCORPORATED 

Frank  H.  Binder 
Program  Manager 
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